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A simple example of a non-equilibrium system for which fluctuations are important is a system 
. of particles which diffuse and may annihilate in pairs on contact. The renormalization group can 

• be used to calculate the time dependence of the density of particles, and provides both an exact 

' value for the exponent governing the decay of particles and an e-expansion for the amplitude of 

this power law. When the diffusion is anomalous, as when the particles perform Levy flights, the 
critical dimension depends continuously on the control parameter for the Levy distribution. The 
e-expansion can then become an expansion in a small parameter. We present a renormalization 
group calculation and compare these results with those of a simulation. 
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Many different approaches have been used to study the dynamics of systems far from equilibrium [j], 0, Q ■ These 
include exact solutions derived by mapping the problem to a quantum spin system Q or by studying various particle 
distribution functions Q, renormalization group methods and simulations. Many of these systems can be 
assigned to a small number of universality classes, based on the time dependence of a few measurable quantities. 
\ Renormalization group methods are particularly useful in characterizing this universal behavior, as they make it 
possible to examine the action describing the behavior of the system and determine which scaling variables are 
relevant. Systems of particles which can diffuse and undergo reactions are straightforward examples, which can be 
studied far from equilibrium. One of the simplest examples of a reaction-diffusion system is the pair-annihilation 
reaction, in which members of a single species of particle, denoted by A, react at some rate A to form an inert 
product. This reaction is written A + A — > 0. This reaction, and some related ones, have been studied for some 
time. For some examples, see the work of Smoluchowski 0, who studied the coagulation of colloidal particles unsing 
O \ a mean- field approach, and Ovchinnikov and Zeldovich [8|] , who examined the effects of fluctuations on the reaction 
A + B — > 0. The solution to the (mean-field) rate equation 

^C) , for the density of particles is the correct result for spatial dimension d > 2, but for lower dimensions, fluctuations 
become important. If the number density of particles scales as n ~ At~ a for large t, then a — 1 for d > 2. A 
renormalization group study by Lee ||| produced the exact result a = d/2 for d < 2, and also yielded the amplitude 
m : a as an expansion in e = 2 — d. The agreement between this amplitude and the exact result [lOj for a specific model 
in d = 1 was poor, as might be expected as the expansion parameter e = 1 is large here. However, the e-expansion 
does provide a systematic picture of the scaling behaviour of this process. 

This paper extends the renormalization group calculation of [fj to the case of anomalous diffusion, which is modelled 
by a long-range hopping process in which the distance a particle travels in each time step is chosen from a Levy 
■ distribution. These distributions form a family that share with the gaussian the property that they are "stable," 
in the sense that the probability distribution of the sum of two numbers chosen from a Levy distribution is the 
same distribution, up to a trivial rescaling. The probability distribution for each of these distributions has a Fourier 
^ ' transform 
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P(k) = e-° AK ~ , (0<cr<2), (2) 



where a is a parameter that controls the shape of the distribution and Da scales the distribution. For a < 2, the 



real-space distributions have power-law tails, P(r) ~ r -(d+<x) f or i ar g e r> These distributions appear in a number of 
physical context s ITlJ . including diffusion in disordered media [lj and motion of particles in turbulent flow, in both 
experiments fl3l Il4j . and theoretical calculations 0] . Hinrichsen and Howard previously simulated the process 
studied here, and determined the exponent a, but did not perform a renormalization group calculation and did not 
measure the amplitude. 
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FIG. 1: The diagrams contributing to the renormalized annihilation rate. Time flows from right to left. 



The calculation of the density of particles uses reasonably standard renormalization group techniques, and will only 
be summarized here. The density of particles can be found by solving the Langevin equation 



dt 



{D N V 2 + D A V ff )^(x, t) - A(/> 2 (x, t) + </>(x, t)C(x, t) 



(3) 



with an initial condition set by the initial density of particles hq. The anomalous diffusion term, ZJ^V^, is defined 
by its action in Fourier space, 



D A V e = -Da k e 



and the noise term £ has correlations of the form 



(C(x, t)({x',t')) = -2X5 d {x - x')5{t - t'). 



(4) 



(5) 



The average density is then given by an average of solutions to Eq. J3J over noise histories, n = (<p). The normal 
diffusion term Dpf~S/ 2 appears if the distribution of hops has any component proportional to k 2 in its Fourier transform. 
This term will be dropped, as it is less relevant than the anomalous term and flows to zero under renormalization. 

The form of the noise term and the noise-noise correlation function given by Eq. JSJ) are not determined by equilib- 
rium physics, but can be derived from the master equation describing the microscopic behavior of the system, using 
the procedure developed by Doi 01 an d Peliti 0] . Their procedure produces an effective field theory describing the 
behavior of the system, with an action given by 
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dt 



D A V a )(l) + 2\(fxt) 2 +\cj) 2 (j) 2 



} - W(0) 



(6) 



The field (j> is a response field, which plays the same role as that in the Martin-Siggia-Rose approach ^^|. This field 
theory can then be used to derive a Langevin equation by integrating out the response field, applying the Martin- 
Siggia-Rose approach in reverse. For examples of this derivation applied to similar problems, see @ and |2Cj . 
Either the field theory or the Langevin equation can be used to develop a renormalized perturbation expansion for 
the density of particles. 

The sum of tree diagrams is equivalent to the solution of Eq. © without the noise term. Divergences appear as 
diagrams with loops are included, and these can be handled by renormalizing the annihilation rate A. Power counting 
shows that the critical dimension is d c — it, so that the expansion will be in e = a — d. The diagrams in Fig. n§i ve 
the full renormalization of the annihilation rate, and represent the sum 
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This sum can be done to all orders after a Laplace transformation A/?(k, s) = L°° dte~ st \r :(k, <), to give 

A 



X R (k = 0,s) = 



l + CT(|)a 



-e/c 



(7) 
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with 



C 



1 2 2 - d /° Y(i) 
a {An) d / 2 F(f) 



(9) 



and so, after a renormalization point s = K a is chosen, the exact flow function for the renormalized annihilation rate 
9r = k.~ c \r(0,k) is 



13 = n 



dg 
d 



- = -eg R + eCT{-)g R . 



(10) 



This has a non-trivial fixed point, which is stable for d < a, at g* R = (CT(-)) -1 . 

Imposing the condition that the density be independent of the (arbitrary) normalization point k, and using dimen- 
sional analysis, a renormalization group equation can be written as 



d d d 

d{D A t) dg R dn 



n R = 0. 



The solution to this, by the method of characteristics, is 

n R (D A t,n Q ,g R ,K) = {k° D A t)- d ^n R ( K -° ,(k° D A t) d '° ,g R , K ) 



(12) 



where g R is the running coupling, which goes to g* R as t — > oo. In a diagrammatic expansion, each loop brings in a 
higher power of the renormalized coupling, so an expansion in the number of loops is an expansion in g R , which is 
small near the critical dimension d c = a. The first approximation to the right hand side can be found by summing 
all tree diagrams generated by expanding either the action (Eq. ©) or the Langevin equation (Eq. J3J)). The next 
term, including all diagrams with loops, is calculated by writing an integral equation for the density, as done in [jj. 
The leading contribution to the density at long times is then given by 



n(t) = A(D A ty 
The exponent is exact, and, to order e°, the amplitude is 



A = 



i r(f)2 1 - d / a 
( 47r )d/ 2 r (|) - 



(13) 



(14) 



This expression is obtained by expanding factors which are singular as e — > in e, and leaving the remainder as an 
expression in d. 

Simulations of a microscopic model in one dimension have been performed and the results compared to the prediction 
of the renormalization group calculation. The convergence of the densities to the predicted power law is clearly shown 
in Fig. H 

The simulations are of a one-dimensional lattice, where only one particle may occupy each site. This differs from 
the renormalization group calculation, in which multiple occupancy is allowed, but the annihilation rate flows to a 
fixed point and thus the bare annihilation rate does not appear in the final answer, so the results should be similar. At 
the beginning of the simulation, the lattice contains L = 10 7 sites, with every site occupied. Whenever the number of 
particles fell below 1000, the system was doubled by appending an exact copy of the current configuration of particles. 
While this does make the system momentarily periodic, the two halves subsequently evolve differently. This allowed 
the simulation to continue to large times without large statistical fluctuations in the density. For a = 1.05, the final 
system size was L = 2 9 x 10 7 . The distribution of jumps was chosen to follow a Levy stable law, using the method 
given in |2l|. Since this method is for a continuous distribution, which is then made into a discrete distribution by 
rounding, the values generated were multiplied by a numerically determined factor to produce the correct low-fc value 
of the distribution. Earlier simulations |l6U20j used a pure power-law form for the distribution of jumps, which does 
flow to the desired distribution in the long-time limit. The distribution used here matches the Levy distribution much 
more closely after a single step, and makes it easier to determine the anomalous diffusion constant D A . 

For some values of a, the anomalous diffusion constant D A in the distribution of hop lengths, Eq. @, was varied. 
In the field theory or the Langevin equation, this change to the kinetics results in a change in the coefficients of powers 
of k, which come from an expansion of the hop length distribution. These terms, many of which are irrelevant in the 
renormalization group sense and are not written in the Langevin equation, can be seen to have a significant impact on 
the crossover behavior, and it may be useful to vary these parameters in future simulations, to examine the crossover. 
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FIG. 2: The density of particles divided by its asymptotic (large time) power law, for several values of a. For a close to 1 
and Da = 1 (open symbols), the crossover becomes very slow, but the trend is towards the same value as for simulations done 
with a larger value of Da- The filled symbols give the density for Da = 3.5. Also shown (solid line) is the crossover function 
of Eq. 115H for a — 1.1, with the parameters extrapolated to their t — > oo values. 



• simulation 
-|-exact result 
e-expansion 




FIG. 3: The amplitude of the power law decay of the density of particles determined in simulation (circles), compared with 
the renormalization group prediction of Eq. C1J. Also shown is the exact result [T(i| | for the normal diffusion case (cross). 



As can be seen in Fig- EI the time taken to reach the asymptotic form of the density depends strongly on Da- For 
a close to 1 and Da = 1, the crossover is very long, but the amplitude can be extracted from a fit to the form 

n{t) A(l-Bt-+). (15) 



To obtain the long-time value of the amplitude, the fit was done over many ranges with differing starting times, and 
extrapolated to t — > oo. For larger values of Da, this fit is not necessary, as the density reaches its asymptotic value 
quickly, and several decades of scaling can be seen in the data. 

The prediction for the amplitude is compared to simulation results in Fig. El The agreement between the e-expansion 
for the amplitude and the simulation result becomes quite good for small e, as expected for this asymptotic power 
series expansion. As well as being interesting as a model for physical processes with anomalous diffusion, the Levy 
flights used here allow this regime, where the expansion parameter is small, to be explored in a simulation. 
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